Surface solitons in two-dimensional quadratic photonic lattices 
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We study two-color surface solitons in two-dimensional photonic lattices with 
quadratic nonlinear response. We demonstrate that such parametrically cou- 
pled optical localized modes can exist in the corners or at the edges of a square 
photonic lattice, and we analyze the impact of the phase mismatch on their 
properties, stability, and the threshold power for their generation. 
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Two-dimensional surface solitons have been recently 
predicted to exist as novel types of discrete solitons lo- 
calized in the corners or at the edges of two-dimensional 
photonic latticesii^ii^. These theoretical predictions 
were followed by the experimental observation of two- 
dimensional surface solitons in optically-induced pho- 
tonic lattices^ and waveguide arrays laser-written in 
fused silica^!. Importantly, these two-dimensional surface 
solitons demonstrate novel features in comparison with 
their counterparts in truncated one-dimensional waveg- 
uide arrays. In particular, in a sharp contrast to one- 
dimensional surface solitons, the threshold power of two- 
dimensional surface solitons is lower at the surface than 
in a bulk making the mode excitation easier^. 

Surface solitons are usually considered for cubic or sat- 
urable nonlinear media. However, multicolour discrete 
solitons in quadratically nonlinear lattices have been 
studied theoretically in both one- and two-dimensional 
lattices^iii^i^ irrespective to the surface localization ef- 
fects. Only Siviloglou et &\}^ studied discrete quadratic 
surface solitons experimentally in periodically poled 
lithium niobate waveguide arrays, and they employed a 
discrete model with decoupled waveguides at the second 
harmonics to model some of the effects observed experi- 
mentally. 

More elaborated theory of one-dimensional surface 
solitons in truncated quadratically nonlinear photonic 
lattices, the so-called two-color surface lattice solitons, 
has been developed recently by Xu and Kivshar— who 
analyzed the impact of the phase mismatch on the ex- 
istence and stability of nonlinear parametrically cou- 
pled surface modes, and also found novel classes of one- 
dimensional two-color twisted surface solitons which are 
stable in a large domain of their existence. 

In this Letter, we extend the analysis of two-color sur- 
face solitons to the case of two-dimensional photonic lat- 
tices. We study, for the first time to our knowledge, two- 
color surface solitons in two-dimensional square photonic 
lattices with quadratic nonlinear response. We analyze 
the effect of mismatch on the existence, stability, and 
generation of surface solitons located in the corners or at 
the edges of the nonlinear lattice. 



We consider the propagation of light in a two- 
dimensional photonic lattice of a finite extent imprinted 
in a quadratic nonlinear medium, which involves the in- 
teraction between the fundamental frequency (FF) and 
second-harmonic (SH) waves. Light propagation is de- 
scribed by the following coupled nonlinear discrete equa- 
tions^iii 
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27"«,m'^n,m exp(-|-i/3z) = 0, 
exp(-i/3z) = 0, (1) 



where it„^m and Vn,m EW'e the normalized amplitudes of 
the FF and SH waves, respectively, Cu and are the 
coupling coefhcients, 7 characterizes the second-order 
nonlinearity, and f3 is the effective mismatch between two 
harmonics. The second-order difference operator A2 is 
defined as 



A2U„ - 
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Equations ([T]) conserve the total power, 

P^P^ + P,^ J2iWn,m\^ + 2|w„,™n. (2) 
n,m 

We look for stationary two-mode solutions of Eq. ([T|) 
in the form, Un,m{z) = Un,m'sx.'9{i^z) and Vn^rn{z) = 
Vn^m exp(2zAz — ifiz), and obtain the nonlinear algebraic 
equations for the (real) mode amplitudes, Un,m and Vn,m, 



= 0, 
= 0, 



(3) 



where for the (NxM) lattice we have n = 0, 1, . . . , TV, 
m = 0, 1, . . . , Af and Un,m = and Vn,m = if either 
n < 0, or TO < 0, or n > TV, or m > M. 

Equations ([3]) possess symmetry properties. For in- 
stance, the transformation 7 —7, Vn.m —Vn^m leave 
Eq. ([3]) invariant. 

In the anti-continuum limit, i.e. when the couplings 
in the lattice vanish, Eqs. ([3]) imply Vn,m = -^/27 and 
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Fig. 1. (color online) Examples of the corner (left) 
and edge (right) localized modes in two-dimensional 
quadratic lattice. Shown are the amplitudes of the FF 
components for Cu — I, Cy — 0.5, (3 — 0, and A = 5. 

L/2 „ = (A/27)(2A - /3)/7. Thus, if A > 0, then 2A > (3. 
If A < 0, then 2A < /3. Also, localized modes should 
exist outside the hnear spectrum band, |A| > 4C„. Thus, 
in the propagation constant-mismatch space, the region 
where localized modes exist is bounded by A > iCu, A > 
(3/2, or A < -4C„, A < /3/2. Also in this hmit, the 
propagation constant is proportional to rather than 
P, as in the case of the cubic nonlinearity. 

We consider a square lattice of a finite extent with 
iV = M = 21, and look for two-dimensional localized 
modes in the corner (1,1) and at the edge {N/2, 1) of 
the lattice. We use two sets of the coupling parame- 
ters: (i) strong couphng, C„ = 1 and Cy = 0.5, this is 
similar to the parametric processes in bulk media, and 
(b) no coupling for the second harmonic, C„ = 1 and 
Cy — 0, similar to the fabricated structures employed in 
the recent experiments^^. For given values of the cou- 
pling parameters (C„,C„), nonlinear parameter 7, and 
mismatch (3, we use a standard numerical procedure of 
continuation from the anti-continuum limit. Examples of 
the nonlinear corner and edge nonlinear localized states 
are shown in Fig. [1] where the fields propagating along 
the waveguides are presented as a superposition of the 
waveguide modes, U{x, y) = U„,m4'ix - n,y - m), 

where (j){x, y) is the (single mode) guide centered on site 
(n,m), and similarly for V{x,y). In Figs. 1 and 4, we 
use 4>{x, y) = exp(-x^ - y^)/a'^, with — 0.05. 

We focus on the analysis of the power dependencies 
characterizing the families of two-dimensional localized 
modes. The typical dependencies of the total power ^ 
on the propagation constant A are shown in Fig. [21 for 
several values of the mismatch parameter (3 marked on 
the plots. The power curve has a typical minimum cor- 
responding to the minimum (threshold) power for the 
localized mode to exist. In the region where the slope of 
the power P becomes negative, the localized modes are 
unstable, as observed for other types of surface modes 
and also verified here by numerical analysis. 

Figure [3] shows the dependence of the minimum power 
Pmin required to create a stable localized mode vs. mis- 
match P, either in the corner or at the edge of the lattice. 
In general, the results for the corner and edge modes are 
rather similar. Clearly, the value of the mismatch de- 
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Fig. 2. (color online) Total power of the two-dimensional 
modes localized in the corner (left) and at the edge (right) 
of the lattice for Cu = 1, Cy = 0.5, and for three different 
values of the mismatch parameter (3. 

termines the threshold power needed to generate a local- 
ized state, with a minimum value observed for positive 
f3 larger than (3 = 5. The minimum power depends also 
upon the value of the coupling constants. However, the 
main result is that for virtually all values of the mismatch 
parameter the corner mode requires less power than 
the edge mode; this is somewhat similar to the case of 
the cubic nonlinear latticeii^ii^. 




-10 -5 5 10 -10 -5 5 10 



yS IS 

Fig. 3. (color online) Minimum power vs. mismatch for 
the corner (circle) and edge (square) modes, for (a) C„ = 
1, Cy = 0.5 and (b) C„ = 1, C„ = 0. For virtually 
all values of the mismatch parameter, the corner mode 
requires less power than the edge mode. 

Using the same numerical approach and starting from 
the anti-continuum limit, we find several other families 
of two-dimensional localized modes located in a close 
vicinity of the lattice corners and edges. These modes 
provide a two-dimensional generalization of the surface 
modes known for one-dimensional lattices placed at dif- 
ferent distances from the edge, and corresponding to a 
crossover between the surface and bulk discrete solitons 
as discussed earlier^. In addition, we find novel classes 
of the so-called two-dimensional twisted modes, an exam- 
ple of one of such modes is shown in Fig. [4]for C„ = 1, 
Cy = 0.5, and A = 5. 

To study the generation of these two-color surface 
modes, we launch a tight beam at one site, Mn,m(0) = 



3 



Uo5n^„gSm,mo and Vn,m{0) = VoS„^n„6m,mo. either in the 
corner or at the edge of the lattice. We let the system 
evolve from 2; = up to z = 2,„ax, and trace the shape 
of the beam measuring the partial powers P„ and Py for 
both FF and SH components, respectively. In general, 
the results are qualitatively similar for the two initial 
beam positions (in the corner and at the edge) . The only 
difference is in the actual value of the field amplitude at 
which the beam self-trapping occurs. On the other hand, 
significant differences are observed between two sets of 
the coupling parameters and also whether both fields or 
only the FF field is initially present. 




Fig. 4. (color online) Example of a two-dimensional 
twisted mode at the lattice corner. Shown are the am- 
plitude of the FF (left) and SH (right) fields for C„ = 1, 
C„ = 0.5, A = 5, and /3 = 0. 




Fig. 5. (color online) Excitation of two-color surface 
modes at the lattice corner. Shown is the dynamics of 
FF (solid) and SH (dashed) components for (a,b) Cu — I, 
Cy ^ 0.5; (3 = with (a) Uq = Vq = 5 and (b) Uo = 20, 
Vo = 0; and (c,d) C„ = 1, C^, = 0, with (c) Un = 5, 
VJ, = 0, /3 = 0, and Uq = 20, Vq = 0, /S = -5. 

In general, when both the fields are initially excited 
and the power is strong enough, a localized surface mode 
is formed, whose amplitude oscillates at a relatively fast 
rate [see Fig. (SJa)]. Presence of a mismatch tends to de- 
stroy the localization, but it is restored if the initial power 
is high enough. When the SH field is initially absent, we 
observe that it is harder to excite a localized mode [see 



Fig. [3b)]. When it is achieved, the mode always seem 
to oscillate between zero and a maximum amplitude at 
a fast rate. When the SH coupling is absent (Ct, = 0) 
and the harmonics are matched (/3 = 0), it can be proven 
that the resulting equation for the FF field is transformed 
into a discrete equation with a nonlinear loss term: 
i{d/dz)un,m + C„A2U„,m -l- 2ij'^u'^ „^ ul „^{s)ds = 0. 
This induces a quick power loss in the FF field and its 
rapid transfer to the SH field [see Fig. [^c)]. The ad- 
dition of mismatch and a high enough initial FF field, 
restores self-trapping for both the fields, where both the 
fields oscillate at a fast rate [see Fig.lSJd)]. 

In conclusion, we have studied localization of light 
in two-dimensional quadratically nonlinear photonic lat- 
tices and determined the conditions for the existence of 
two-color surface states localized in the corners or at the 
edges of the lattice. We have analyzed the impact of 
the phase mismatch on the properties and stability of 
two-dimensional localized modes, as well as the thresh- 
old power for their generation. 
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